Critical Temperature of a Trapped Bose Gas: Mean-Field Theory and Fluctuations 
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We investigate the possibilities of distinguishing the mean-field and fluctuation effects on the 
critical temperature of a trapped Bose gas with repulsive interatomic interactions. Since in a direct 
measurement of the critical temperature as a function of the number of trapped atoms these effects 
are small compared to the ideal gas results, we propose to observe Bose-Einstein condensation 
by adiabatically ramping down the trapping frequency. Moreover, analyzing this adiabatic cooling 
scheme, we show that fluctuation effects can lead to the formation of a Bose condensate at frequencies 
which are much larger than those predicted by the mean-field theory. 

PACS number(s): 03.75.Fi, 67.40.-w, 32.80. Pj, 42.50.Vk 



The experimental observation of Bose-Einstein con- 
densation in atomic gases of *^Rb ^"^Na [||, and ^Li 
confined in a magnetic trap, has enormously boosted 
the theoretical and experimental research of such quan- 
tum gases. After studying the elementary excitations of 
the condensate [Q,D , the JILA and MIT groups have also 
started to determine more accurately the critical temper- 
ature of the e clS 3tS cl function of the number of atoms in 
the trap, and in particular the deviation of this temper- 
ature from the ideal case . Measuring these deviations 
will give important information on the effect of interac- 
tions on the thermodynamic properties of the gas and 
may also signal a break-down of the Bogoliubov theory 
that has been so succesful in explaining the collective ex- 
citation spectrum |]7|-|l0|]. It is therefore of considerable 
theoretical interest. 

For a non-interacting gas in an isotropic external har- 
monic potential V{r) = ma;^r^/2, it is well known that 
if the number of particles N ^ 100, Bose-Einstein con- 
densation occurs at a temperature 
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However, interactions between the atoms influence the 
critical behavior of the alkali gases of interest and shift 
the critical temperature. Roughly speaking there are two 
effects to consider. First, the inhomogeneity of the gas 
leads (above the transition) to a spatially-varying mean- 
field potential of 2n{r)T^^ , where n(r) is the local den- 
sity of the gas, T^^ ~ Anah^ /m is the two-body tran- 
sition matrix at zero energy and momenta, and a is the 
s-wave scattering length. Depending on the sign of the 
scattering length the atoms are therefore either repelled 
from (a > 0) or attracted to (a < 0) the center of the 
trap and it requires more or less atoms respectively, to 
acquire the necessary central density for Bose-Einstein 
condensation. Using a local-density approximation this 
effect was recently studied by Giorgini, Pitaevskii and 



Stringari They obtained a mean-field shift in the 

critical temperature of 



Tn 
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where I — ^fij muj is the extent of the harmonic oscil- 
lator ground state. Secondly, Bijlsma and Stoof have 
found by means of a renormalization group calculation 
that for gases with a positive scattering length the nec- 
essary central density for Bose-Einstein condensation is 
reduced considerably due to fluctuation effects |l2| . This 
means that the degeneracy parameter n(0)A^ (with 
At,, = {I'n'h^ ImkBTfl'^ the thermal de Broghe wave- 
length) at the critical temperature is not C(3/2) — 2.612, 
as assumed by the mean-field theory, but instead a mono- 
tonically decreasing function of the quantum parameter 
a/Kth- It is therefore always smaller than the ideal gas 
value C(3/2). 

Although the latter effect clearly increases the criti- 
cal temperature, it is at present unclear to what amount 
this increase cancels the mean-field shift in the critical 
temperature. To investigate this issue, we need quan- 
titative information on the density profile of a trapped 
Bose gas with repulsive interactions, which we obtain by 
performing numerical calculations in the local-density ap- 
proximation. In this approximation the system is consid- 
ered to be locally homogeneous, and the inhomogeneity 
induced by the trapping potential V{v) and the mean- 
field interaction 2n(r)T^^ is determined selfconsistently 
by the spatially-dependent effective chemical potential 
^'(r) = ^ — y(r) — 2n(r)T^^, where ^ is the actual chem- 
ical potential of the gas. 

For our purposes, which will become clear shortly, we 
also need to calculate the total number of particles N in 
the gas, together with the entropy S. According to sta- 
tistical physics, these quantities can be calculated from 
the grand-canonical potential $7 of the system, i.e. from 
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N 



and S = 



T.V 



df 



Above the critical temperature the (mean-field) grand- 
canonical potential of an interacting Bose gas is given 
by 

n^keTY, - cxp[ - Pie, - fi)]) 

^ j dv n'{v)T^^ , (3) 

where denote the one-particle energy levels in the ef- 
fective potential y(r) -I- 2n(r)T^^ Implementing the 
local-density approximation by putting = ?i^k^/2m -|- 
V{r) + 2n(r)T2'B and using Y.j = I dr J dk/(27r)3, we 
can perform the gaussian integral over the momenta Kk 
to find 

- J dr 35/2 [2 (r)] 

for the first term on the right-hand side of Eq. . Here 
we introduced the Bose function 
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and the local fugacity z{r) = exp [(3^'{r)]. 

Calculating now the derivative of with respect to /x, 
one obtains for the total number of particles 



N = J dr n{r) = J dr 



93/2[z{r)], 



(5) 



where we used 53/2 (-2) = z[dgc,^2iz)/dz]. Performing a 
similar differentiation with respect to the temperature, 
the entropy of the system is found to be 



dr 



;95/2[z{r)] - nir) ln[z(r)] 



(6) 



Note that in both differentiations we used the fact that 
z{r) has an implicit dependence on the chemical potential 
and the temperature through the density n{r). The con- 
tributions to the particle number and the entropy from 
this dependence exactly cancel, however, the contribu- 
tions coming from the second term in the right-hand side 
ofEq. (|). 

At this point we are ready to discuss the way in which 
we numerically determine the critical temperature of the 
gas as a function of the number of particles in the trap 
and, in particular, how we include the fluctuation effects 
on this critical temperature. To include these fluctua- 
tions in a first approximation, we calculate at a given 
temperature the value of the quantum parameter a/Att 
and extract from the appropriate value of the critical 
degeneracy parameter the value of the central den- 
sity n(0) at which Bose-Einstein condensation first oc- 
curs. Knowing the critical central density we can also im- 
mediately determine the chemical potential at criticality 




from n(0)A3^ = 33/2 [exp - 2n(0)r2S)}] xhe cor- 
responding number of particles is then self-consistently 
calculated from Eq. ^ . To compare with the mean-field 
results, we also calculate the number of particles in the 
system for a central density of C(3/2)/Aji_^ Both results, 
together with Eq. (^) are plotted in Fig. |l|. 



0.00 



-0.02 



P= -0.04 



-0.06 



-0.08 

0.00 0.02 0.04 0.06 

(a/DN"" 

FIG. 1. Relative change in the critical temperature for an 
interacting gas in a harmonic oscillator trap. Curve 1 gives 
the mean-field result and curve 2 includes also the effect of 
fluctuations. The dashed line represents Eq. (2). 

For different values of the scattering length a and the 
trapping frequency w, we determined the critical temper- 
ature as a function of the number of particles and found 
that for both calculations the relative shift in the criti- 
cal temperature is a universal function of {a/l)N^/^ as 
indicated in Fig. From this figure wc see that Eq. (|^) 
gives the correct lowest order mean-field shift of the tem- 
perature, but also that higher order corrections become 
important relatively soon. Moreover, the fluctuation ef- 
fects are roughly speaking of the same order of magni- 
tude under the experimental conditions of interest, al- 
though they are clearly unable to overcome the decrease 
of the critical temperature due to the mean-field poten- 
tial. However, at this point we have to remind ourselves 
that the present calculation underestimates the infiuence 
of fluctuations because it does not incorporate the effect 
that fluctuations have on the density proflle in the trap. 
In principle, this could also be included by performing 
what one might call a 'local renormalization group' cal- 
culation, i.e. apply the renormalization group theory of 
Ref. at each point in space. 

Since the critical temperature of an interacting Bose 
gas in an absolute sense does not deviate much from the 
ideal gas value, we now propose a possible experiment 
that could make the effects of the interatomic interac- 
tions visible without having to deal with a large ideal gas 
'background'. Suppose we have a gas of N particles in 
an isotropic trap with frequency uj and at a temperature 



2 



T above the critical temperature. The gas also has a cer- 
tain entropy given by Eq. . We now ask ourselves what 
would happen if the system is cooled adiabatically to a 
lower temperature by slowly lowering the trapping fre- 
quency to a value lo' . Since during this adiabatic process 
the entropy and the number of particles in the gas re- 
main constant, we know the answer to this question if we 
determine the final chemical potential and temperature 
combination {^i' , T') such that iV(/x', T', lo') = 7V(/i, T, w) 
and S{ii',T',Lu') ^ S{fi,T,uj). 

The reason why we propose to adiabatically cool the 
gas in this manner is that for the ideal Bose gas, the 
only way to keep both the entropy as well as the par- 
ticle number constant during the change u u' in 
the trapping frequency, is to change the temperature 
from T to {lo' /lo)T and the chemical potential from ji 
to (lo'/lo))!. This leaves all the occupation numbers of 
the harmonic oscillator states unchanged. More impor- 
tantly, it leaves the fugacity and therefore the degeneracy 
parameter n(0)Aj^ in the center of the trap unchanged. 
So for the ideal gas, adiabatic cooling does not bring 
us closer to criticality. Therefore, if we can reach Bose- 
Einstein condensation in an interacting gas in this way, 
it must be entirely due to the effect of interactions. 

In Fig. H, the value of the degeneracy parameter in the 
center of the trap is plotted for five different adiabatic 
processes (solid lines 1 to 5) during an adiabatic ramp 
down of the trap frequency from a starting frequency of 
373 Hz. All five adiabatic curves correspond to a gas 
containing N = 20,000 ®'^Rb particles with an s-wave 
scattering length of 109ao. 
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FIG. 2. Adiabatic trajectories (solid lines) of a gas of 
20, 000 *^Rb particles when the trapping frequency is ramped 
down. The initial temperatures are 1) 442.21 nK, 2) 445.85 
nK, 3) 449.43 nK, 4) 454.72 nK, and 5) 464.16 nK. The dashed 
line indicates the value for the central degeneracy parameter 
at which Bose-Einstein condensation first occurs if fiuctuation 
effects are taken into account. The dot-dashed line denotes 
the mean-field critical degeneracy parameter. 



The initial temperature of the gas increases from the 
upper to the lower curves and consequently the chemi- 
cal potential increases from the lower to the upper curve, 
since the total number of particles is the same. Physi- 
cally, for a nonideal gas the extra term 2n{r)T^^ /ksT 
in the exponent of z(r) causes a deviation from a hor- 
izontal line, because the density does not scale linearly 
with T but approximately with T^^^. Thus keeping N 
and S constant when lo is varied, we now have to change 
the local fugacity and as a result the degeneracy param- 
eter in the center of the trap also changes. Furthermore, 
to achieve Bose-Einstein condensation in a gas with a 
positive scattering length, it is clear that the frequency 
must be ramped down because this lowers the density of 
the gas and reduces the effect of the repulsive interac- 
tion 2n{r)T^^ . As a result, the density in the center of 
the trap reduces less than in the ideal gas case and the 
degeneracy parameter n(0)Aj^ increases. 

In Fig. g we also plotted the value of the critical de- 
generacy parameter following from the renormalization 
group theory of Ref. and corresponding to the in- 
stantaneous temperature during the adiabatic cooling of 
the gas. It should be noted that in principle this curve is 
different for each adiabatic trajectory. However, for the 
five trajectories presented here the differences cannot be 
distinguished on this scale. Another aspect that can be 
seen from Fig. ^ is that only the upper three curves reach 
the conditions for Bose-Einstein condensation, whereas 
curves 4 and 5 remain uncondensed all the way down to 
u! = which is equivalent to T = 0. Indeed, a detailed 
analysis of the chemical potential and the central density 
as a function of the instantaneous temperature T shows 
that for curves 1 to 3 /it oc T" and n{0) (x for the com- 
plete cij-interval considered. Moreover, the powers obey 
a < /? < 3/2 which implies that although initially the 
gas is not condensed and fi < 2n(0)T^^, there always 
will be some lower and nonzero temperature such that 
/i = 2n(0)T2^ and the gas fulfills the mean-field critical 
conditions. Consequently it will have already achieved 
the critical conditions including fluctuation effects at a 
higher temperature. 

However, for curves 4 and 5 the chemical potential and 
central density exhibit a different behavior. This is due 
to the fact that during the cooling process the chemical 
potential becomes negative in order to satisfy conserva- 
tion of both the number of particles A'' as well as the 
entropy S. Therefore, the system cannot become Bose 
condensed unless ^ t and 2n{0)T'^^ i for T 0. 
Looking again in more detail to the power behaviour of 
/J, and n(0) as a function of T, wc find that n{0) cx 
with P almost equal to 3/2, and that /i oc —T for T ^ 0. 
So both fi and 2n(0)T^^ indeed approach zero, but the 
limit of the ratio fj,/kBT will now decide whether or not 
the system achieves Bose-Einstein condensation at zero 
temperature and frequency. 

Fig. 1^ shows a plot of the trapping frequencies at which 
Bose-Einstein condensation first occurs if only mean-field 
effects are taken into account (curve 1) and if also fluctu- 
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ation effects are included (curve 2). Again all adiabatic 
processes are for a gas containing 20, 000 particles and 
an initial trapping frequency of 373 Hz. The initial tem- 
perature of the gas is given by the ordinate. From this 
figure it is clear that the corrections due to fluctuations 
have a large influence on the flnal frequency at which 
Bose-Einstein condensation is achieved. For suitable ini- 
tial conditions, the flnal trapping frequency can exceed 
the mean-field prediction by more than a factor of 3. We 
therefore conclude that in principle the fluctuation cor- 
rections are experimentally observable by adiabatically 
cooling the gas. 

400 I ■ . ■ . ■ 1 
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FIG. 3. Trapping frequency needed to obtain the critical 
conditions in the center of the trap, starting from a configura- 
tion of 20,000 particles in a 373 Hz trap and at a temperature 
denoted by the ordinate. Curve 1 is the mean-field result and 
curve 2 includes also the effect of fiuctuations. 

However, it must be noted that while curves 1 and 
2 in Fig. 3 can be separated by as much as a factor of 
three in the frequency axis, they are separated by only 
1% on the temperature axis. This should be compared 
with the JILA group's most recent attempt at absolute 
thermometry, with claimed absolute accuracy no better 
than 5% Nevertheless, measurements of the conden- 
sate fraction near the critical temperature indicate that 
the shot-to-shot reproducibility of the scaled temperature 
is much better, i.e. on the order of 1% [Q. Moreover, 
the scaled temperature varies smoothly with flnal evap- 
orative rf cut, which is a readily adjustable experimen- 
tal parameter. We therefore envision an experiment in 
which samples are prepared at temperatures precisely de- 
termined relative to the empirical onset temperature of 
BEC, and then subjected to adiabatic ramps of the trap- 
ping frequency. 

In summary, we showed that, although fluctuation cor- 
rections decrease the critical density of a trapped gas of 
repulsively interacting bosonic atoms, they do not com- 



pletely cancel the increase in the critical number of par- 
ticles obtained from mean-fleld theory. The effect of the 
interactions, and in particular, of the fluctuation correc- 
tions to the mean-field theory can be observed more di- 
rectly by adiabatically cooling the gas to criticality. Here 
we imply changing adiabatically the frequency of a har- 
monic trap and not the power law of the trapping poten- 
tial, which has recently been shown to change already the 
degeneracy parameter of the ideal Bose gas Q . In addi- 
tion to being more difficult to implement experimentally, 
the latter approach is not so favorable for our purposes 
because it is less sensitive to the effect of interactions. 
We hope that this work might stimulate an experimental 
observation of these interesting phenomena in the near 
future. 
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